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The Majorana phase boundaries in planar 2D hybrid (semiconductor-superconductor) nanowires 
are modified by orbital effects due to off plane magnetic components. We show that Majorana zero 
modes survive sizable vertical field tiltings, uncovering a remarkable phase diagram. Analytical 
expressions of the phase boundaries are given for the strong orbital limit. These phase boundaries 
can be fulfilled with attainable setups, such as an InAs nanowire of 150 nm in transverse width. 
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I. INTRODUCTION 

The physics of 2D electron gases in magnetic fields has 
proved invaluable for the Condensed Matter field with, 
e.g., the celebrated quantum Hall effecti as well as with 
many devices based on quantum wells, wires and dotsi^ 
On the other hand, Majorana zero modes in quasi ID 
wires have recently attracted strong interest, both as 
exotic quasi-particles and as candidates for topological 
quantum computingj^^— In this article we relate 2D-gas 
properties and Majorana physics, showing the remark¬ 
able role of the orbital motions characteristic of 2D sys¬ 
tems in magnetic fields. 

Majorana modes in quasi ID wires are effectively 
chargeless, zero-energy quasiparticles. They arise from 
the splitting, through a phase transition, of bulk elec¬ 
tronic states into pairs of quasiparticles on the wire ends, 
each one being its own antiparticlei^r— An important 
feature of the Majorana mode is that it appears only 
when a critical value of the external magnetic field, a 
phase-transition threshold, has been surpassed. Several 
experiments with hybrid superconductor-semiconductor 
nanowires using tunneling spectroscopy from a normal 
conductor to the nanowire have observed a zero bias peak 
consistent with a Majorana stateJ^^— The observed peak 
height is, however, an order of magnitude lower than the 
quantized value 2e^//i. This discrepancy is not yet well 
understood, as it might be due to effects ranging from fi¬ 
nite temperature, experimental and tunneling resolutions 
to other low-energy subgap states and possible inelastic 
and renormalization processesi^^iSi 

In practice, distinguishing zero-bias peaks due to Ma¬ 
jorana modes from other potential sources relies on the 
detailed knowledge of the phase diagram in each partic¬ 
ular physical realization. Therefore, it is highly relevant 
knowing how Majorana physics is affected by the extra 
dimension in 2D, with respect to ID. This question has 
been addressed with quasi-lD multiband wiresj^^— mod¬ 
elling transverse modes as mutually coupled ID wires. 
This leads to essentially the ID physics with only one 
Majorana mode at each nanowire end when an odd num¬ 
ber of transverse modes are above their critical magnetic 
field. However, the role of the magnetic orbital mo¬ 



FIG. 1. Schematic of a 2D planar nanowire showing the axis 
definitions. A magnetic field in a tilted direction is included. 
The density distribution of Majorana modes on the wire ends 
is qualitatively shown. 

tion has been usually disregarded. Addressing Majorana 
physics in 2D systems with orbital motion is therefore rel¬ 
evant as a way to discard alternative scenarios that have 
been suggested, such as attributing the observations to 
Kondo-like interaction effectsi^ We want to clarify too 
that the above mentioned experimentsi^r— used cylin¬ 
drical nanowires, for which the orbital effects may be 
different from those discussed here. 

In this article we show that in a planar nanowire the or¬ 
bital motion strongly affects the phase transition bound¬ 
aries. We show how the Majorana phase survives sizable 
vertical field tiltings (Fig. [Ij, even reaching the purely 
perpendicular orientation in some cases. This is not an 
intuitive result since the electronic orbital motion might 
lead to a gap closing, allowing edge propagating solu¬ 
tions, or it might totally change the character of the 
topological states. In this sense we note that the s-wave 
nanowires inside a magnetic field have an associated di¬ 
rectionality with their Majoranas located at the left and 
right edges (Fig.[T|), while a 2D p-wave nanowire has their 
edge states located all along the ID perimeter. 

In the strong orbital limit, the phase transitions occur 
for critical values of the polar angle, following a simple 
analytical law that does not depend on sample details. 
With parallel (x) field orientation the transition law is 
also analytical, while for intermediate regimes the phase 
transitions are obtained numerically. We assess the con¬ 
sistency between the phase diagram and direct calcula¬ 
tions of the Majorana modes in semi-infinite and finite 
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2D wires, emphasizing the importance of covariant grid 
discretizations for the latteri ^°i^^ 

The work is organized as follows. Section II introduces 
the physical system in detail while Sec. Ill contains the 
algorithm used to determine the phase-transition bound¬ 
aries. In Sec. IV we discuss the phase diagrams for 
out-of-plane tiltings of the magnetic field with different 
strengths of the spin-orbit coupling. Sections V and VI 
deal with the actual Majorana solutions in semi-infinite 
and finite 2D Majorana nanowires, respectively. Section 
VII contains the conclusions of the work. 


II. PHYSICAL MODEL 


The combination of s-wave superconductivity, Rashba 
interaction and an external magnetic field is a well known 
source of Majorana Fermionsji We consider electronic 
motion restricted to the x (longitudinal) and y (trans¬ 
verse) directions in presence of these three effects. The 
homogeneous magnetic field points in an arbitrary direc¬ 
tion and the edges are modeled as infinite square well 
potentials in the longitudinal and transverse directions 
(Fig. HI). 

The nanowire physics is described by a Hamiltonian of 
the Bogoliubov-deGennes kind, split in the following way 

'HbcIG ='Hq + 'Hz + Hr + Horb ■ (1) 

The successive contributions to Eq. o are the zero-field 
and superconducting energies. 

-I- V(x, y) - T:, + As ; (2) 

the Zeeman term 
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Hz = Ar (sin 0 cos (jxjx + sin 9 sin (fxTy + cos 9 Oz) \ 
the Rashba coupling term 

'Hr ~ ( Px^y Py^x ) '^z ) 
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and, finally, the magnetic orbital terms 
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In Eqs. (HI)-® we used the following Nambu-spinor con¬ 
vention, relating discrete components with spin (fj-) and 
isospin (trJj.) as 4 'ri, -4'tk)^- 

The contributions in Eq. ([2]) are, in left to right order, 
the kinetic, electrical potential V, chemical potential y 
and superconducting energies. The Pauli operators 
for isospin (particle-hole) are represented by r^^y^z while 
those for spin are (Jx,y,z- The superconductor term repre¬ 
sents an effective mean field approximation to more com¬ 
plicated interactions with a nearby s-wave superconduc¬ 
tor. The Zeeman term, Eq. depends on parameter 


Ar and models the coupling of the spin with a magnetic 
field of arbitrary polar and azimuthal angles (0, (p) = h. 

The Rashba coupling Eq. o is the result of the self¬ 
interaction between the quasiparticle spin with its own 
motion. This interaction is due to the presence of a trans¬ 
verse electric field representing an internal asymmetry in 
the confinement along z that may be either intrinsic or 
externally induced. The first Rashba contribution, de¬ 
pending on Px(^y is called the ID Rashba term while the 
second one, PyC^x, is the Rashba mixing term. 

The joint effects of superconductivity, Zeeman and ID 
Rashba terms give rise to independent Majorana states, 
one from each transverse band like in the ID model. Each 
one of the modes has a different critical magnetic field 


a(c) _ 
— 


{y-enf + A1 


-, 1/2 


, with n = 1,2,, and Cr, 


the transverse mode energies. Adding the Rashba mix¬ 
ing term to this scenario changes the critical fields due to 
the coupling between different transverse bands. It effec¬ 
tively allows only one Majorana zero mode in parameter 
regions where the ID Rashba term would yield an odd 
number of them (even-odd effect) This is further dis¬ 
cussed in the Appendix as a particular analytical limit of 
the most general case presented below. Besides, in ID 
wires the effect of tilting the magnetic field implies the 
additional requirement of the so-called projection rule 
Ar sm9 sincj) < This is due to the indirect gap 

closing of the infinite wire energy bands at ±kc due to 
the tilted field, where kc is a non vanishing wave num¬ 
ber. As discussed in Sec. IV, in 2D nanowires we find 
strong modihcations of the critical magnetic fields, but 
the projection rule still applies. 


In a planar nanowire the perpendicular component 
of the magnetic held induces orbital motions of the 
nanowire quasiparticles. The magnetic orbital terms, Eq. 
([5)1 , describe this motion and their effect on the Majorana 
states is the central point of this article. These terms de¬ 
pend on the magnetic length Iz, dehned as = hc/eBz, 
and they stem from the kinetic and Rashba energies with 
the magnetic substitution Px ^ Px — fiy/^z adding 
the required Pauli matrix for proper particle-hole sym¬ 
metry. In Eq. ® we assumed the Landau gauge centered 
on Pc = 0, although our results are independent on this 
choice as discussed below. 


All parameters of the complete Hamiltonian are con¬ 
stant inside the nanowire, modeled as a perfectly conhn- 
ing box with 'S> Ly. The numerical results of this 
work are presented in characteristic units of the problem 
obtained by taking h, m and the width of the nanowire Ly 
as reference values. That is, our length and energy units 
are, respectively, Ly = Ly and Er = H^/mL^. A spin- 
orbit length Lso is usually defined as Lso = 9?/ma but, 
as explained below, the Hamiltonian orbital terms will in¬ 
troduce an effective transversal confinement. Therefore, 
the comparison between its characteristic length and the 
one of the nanowire real confinement is relevant. We no¬ 
tice that in our convention, the numerical value of a is 
precisely the ratio of transverse and spin-orbit lengths. 
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a/EjjLij = Ly/Lso- 


III. THE MATCHING METHOD 

In topological systems it is in general possible to re¬ 
late the states of the bulk with those at the boundaries, 
a consequence of the bulk-to-edge correspondence prin¬ 
ciple. In our particular case this means that the semi¬ 
infinite Majorana solution 'k will be expressed as a linear 
superposition of the infinite-nanowire eigensolutions 
(i.e., for the same Hamiltonian but disregarding the left 
and right edges), 

k 

where spin and isospin variables are indicated with 
and rjTt respectively. Notice that the infinite-nanowire 
solutions are characterized by a wave number k that, 
accounting for evanescent waves, may be a complex quan¬ 
tity. In a given spin-isospin basis, Xs^{v<y) ^-nd Xs-riVr), 
with Sg. = i and St = ±, the infinite-wire solutions read 

¥'^\x,y,r]^,rjr) = ^ (y) (??a) Xs, (»?r) , 

<S(T Sf 

(7) 

where (y) is a ID 4-component spinor characteristic 
of the infinite-wire solution with wave number k. 

It has been demonstrated that a Majorana phase tran¬ 
sition occurs in a semi-infinite nanowire when the propa¬ 
gating bands for the corresponding infinite nanowire per¬ 
form a gap closing and reopening when increasing the 
magnetic field, at vanishing energy and wave number 
Therefore, to determine the phase boundaries we only 
need to investigate the band structure at k = 0 and 
E = 0. However, a full determination of the band spec¬ 
trum for every set of Hamiltonian parameters by diago- 
nalization is time consuming and computationally inef¬ 
fective. In accordance with this, it has been pointed out 
that in spite of the non locality of the topological states 
a full knowledge of the band spectrum is not necessary in 
general to determine the phase of a topological system, 
but the only relevant regions are those near the Dirac 
cones that appear at the phase transitionsIn our case, 
this implies searching the solutions of 

(s<TSr|/i|s'<Ts'r) = 0 ^ (8) 

s' (jS't 

where h is obtained neglecting all p^j-dependent terms in 
Eq. (dl, 

h= ^^+V{y) - ^ PyCTxTz 


Notice that with Eq. (jS]) we achieved a reduction to an 
effective ID problem. 

We can use the algorithm devised in Refs. HI andto 
study the solutions of Eq. (|S]). The particular parameter 
sets allowing such a solution will signal the gap closing 
of the original 2D Hamiltonian of Eq. o and thus the 
phase transition we are looking for. The algorithm con¬ 
sists in solving Eq. in a ID grid as a linear system, 
assuming vanishing boundary conditions at y = 0 and 
y = Ly. However, due to the homogeneous character of 
this linear system the trivial solution = 0 is al¬ 

ways possible. The algorithm discards the trivial solution 
by means of a matching point ym, where for an arbitrary 
pair of components (so-Sr) = {st) a non vanishing wave 
function is imposed. In addition, continuity of the first 
derivative at the matching point is also imposed for the 
components other than {st), 

= I, (10) 

^ J ^ 

where SEL) j denote grid derivatives using only upper 
(U) or lower (L) y-grid neighbors. 

Equations (nni and m are used at ym in place of 
the Bogoliubov-deGennes ones. In particular, Eq. m 
makes the system no longer homogenous and such that it 
always admits a solution. The algorithm does not ensure 
the first-derivative continuity for the component {st) at 
the matching point. Therefore, this condition is used to 
distinguish the physical from the non-physical solutions 
with the continuity measure 


T = 






( 12 ) 


As mentioned, the E zeroes will signal the desired gap 
closing boundaries of Eq. ©• Further details about the 
algorithm can be found in Ref. [H. 


IV. PHASE DIAGRAMS 

Figure [5)3 shows the phase transition boundaries ob¬ 
tained with the matching method for an Ly = 150 nm 
nanowire with material parameters typical of InAs and 
a magnetic field strength between 0 and 6 T {As = 
0 — 25 Elf). Phase boundaries are signaled by zero values 
of E (red lines). Only E — 0 values represented by red 
lines have physical sense and the color scale is only indi¬ 
cating the measure deviation from zero. As mentioned, 
the energy unit scales as L~^, such that in a 300 nm wire 
25Eu would correspond to 1.5 T. The other panels of 
Fig. [2] correspond to lower ([^ti) and higher (llt,d) values 
of the Rashba coupling strength. We will explicitly calcu¬ 
late the zero modes for particular sets of parameters be¬ 
low. Here, let us anticipate that the orbital terms do not 
change the even-odd effect of multiband nanowiresi^^— 
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FIG. 2. T measure in a color (gray) scale as a function of 
As and polar angle 6. The azimuthal angle remains 0 = 0. 
Zero values of J- (red) signal the phase transition boundaries. 
Phases with a Majorana mode are labeled with an M. From 
top to bottom the panels correspond to a = 0.1 EuLu (a), 
2EuLu (b), ttEuLu (c), AEuLu (d). Dashed lines indi¬ 
cate the analytical limits. We have assumed As = 3 Eu 
and typical InAs parameters g = 15, m* = 0.033. Panel 
(b) corrresponds to an InAs nanowire with a = 30meVnm, 
As = 0.3 meV and Ly = 150 nm in a magnetic field range 
from zero to 6 T. 


The phases in Fig. [5] contain either no Majoranas or 
at most one Majorana mode in regions labeled with an 
M. The main result of this article is that orbital effects 
do not destroy Majoranas into other phases. However, 
they do lead to characteristic phase maps where the Ma¬ 
jorana states survive sizable vertical tiltings (even up to 
0 = 0 in Fig. [5Jl). The change in transition boundaries 
is caused by two reasons. First, by the change in the ef¬ 
fective transversal confinement due the first two terms of 
Eq. ([5]). Note that the first contribution leads to a har¬ 
monic confining. Second, the third term of Eq. ([^ can be 
understood as an effective inhomogeneous magnetic field 
pointing in y direction due to the combination of Rashba 
and orbital effect. 

It is possible to give analytical expressions of the phase 


boundaries in particular limits (see Appendix). For 9 = 
90° and 0 = 0 the critical magnetic fields read 

- en + + A2 , (13) 

where n = 1,2... and e^, are the (transverse) square well 
eigenenergies. This analytical result extends recent find¬ 
ings from other author a^^i^^ who assumed that the contri¬ 
bution in parenthesis in Eq. (US is an effective chemical 
potential from subband n, without specifying its detailed 
a dependence. Analogously, in the strong orbital limit 
As >> As,/r) the critical angles are 

(liTi)) ■ 

where g is the gyromagnetic factor and m* = mime the 
ratio between the electron effective and bare masses (m 
and me respectively). In this limit quasiparticles are con¬ 
fined by the effective A^-dependent harmonic potential 
caused by the first and second terms of Eq. inde¬ 
pendently of the real nanowire transversal boundaries. 
Quite remarkably, Eq. m is independent of the mag¬ 
netic field, Rashba and superconductivity strengths, as 
well as on the specific wire width Ly. In this sense, the 
critical angles are rather universal. The values of Eqs. 
m and (fHl) for the particular parameters used in Fig. 
[^are overprinted as vertical and horizontal dashed lines, 
respectively. As shown, the numerical values match the 
analytical ones in their corresponding limits. In Figs. 
[3r ,b the transitions boundaries do not deviate substan¬ 
tially from the analytical laws in all the plot, while Figs. 
[l:,d show large differences for intermediate values of the 
parameters. 

The structure of the phase transitions in Figure [2^ is 
typical for cases when the kinetic orbital effects already 
dominate around the first transition boundary with in¬ 
creasing Ab. The phase boundaries just bend from a 
vertical to a horizontal line due to the effective change of 
the transversal confinement in the nanowire, from square 
well to a Ag-dependent harmonic confinement [Eq. ([5])]. 
Assuming /i ft! 0, the conditions for this simpler phase 
diagrams {Iz shortest scale) can be written as a triple 
inequality 


(c 


COS 9l 


(c) 




a 


’ EuLf, 



(15) 


Figure [21 d (InAs with Ly = 150 nm) still presents a 
phase diagram qualitatively similar to Fig. although 
the second inequality of Eq. m is not well satisfied (see 
Tab. U in Appendix). On the other hand. Figures ^ 
and d show the modifications of the phase diagram as a 
increases in effective units. As anticipated, the deviations 
can even allow a Majorana state in perpendicular field 
{9 = 0) in Fig. [2Ji. However, so large spin orbit strengths 
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give rise to complicated phase maps that strongly deviate 
from the analytical limits. 

The strong-a effects seen in the lower panels of Fig. 
[2] are caused by the term —ayayfl^ of Eq. ([S]). In¬ 
deed, this term effectively adds a component along y to 
the magnetic held. Therefore, the effective angle 9e is 
such that 6e > 9, thus explaining the downwards shift 
of the lower phase transition boundary in Figs. HJ: and 
d. Figures [2^, b do not change with other azimuthal an¬ 
gles while in the strong-a diagrams (Fig. [2t,d) (j) mod- 
ihes the precise boundary positions for 9 = 90° (ver¬ 
tical lines), but not the horizontal asymptotes and the 
overall qualitative behavior. Note that, as mentioned 
in Sec. m with (p ^ 0 there is an additional require¬ 
ment for the existence of Majorana modes, the projection 
rule Ab sin0sin(() < However, the effective tilt¬ 

ing of the magnetic field towards y caused by the term 
—ayuy/ll does not modify the projection rule because 
this effective tilting is in opposite directions for positive 
and negative y’s while the projection rule applies only to 
homogeneous magnetic fields. 


V. SEMI-INFINITE NANOWIRES 



Explicit zero-energy eigenstates can be obtained in 
semi-infinite and finite nanowires in order to confirm the 
above phase diagrams. We have checked that either one 
or none Majoranas are obtained in the corresponding re¬ 
gions of Eig. [5] In the semi-infinite system exact zero- 
energy solutions can be studied with the complex band 
structure method of Refs. [13 andlslL The calculation is 
carried out solving the boundary condition 

^CfcVl/W,^(y)=0, (16) 

k 

where the allowed complex wave numbers {A:} and 
the 'I'i^s^(j/) functions are obtained with the matching 
method discussed previously in Sec. uni with the only 
difference that k now is not vanishing. 

Equation m can be reworked into 

J2-^k'kCk = 0, (17) 

k 

with the matrix 

Mk'k Jdy (y ). (18) 

Equation inn shows that, when enough wave numbers 
are included, each Majorana state is represented by a 
null-space eigenvector of matrix A4. In Fig. |3^ we can 
see the convergence of the A4 eigenvalues with the cut off 
in wave number for a particular point of Fig. [5]:. Clearly, 
the lower eigenvalue vanishes asymptotically indicating 
that for this point of the phase diagram a Majorana mode 
is present as expected. In Fig. [3 Jd we can see the corre¬ 
sponding density function, confirming the edge character 
of the mode, as also expected for a Majorana. 


FIG. 3. a) Evolution of the lower eigenvalues of matrix M 
when increasing the number of evanescent modes, as given by 
a cutoff |fc|. b) Majorana density function in a semi-infinite 
nanowire for the null eigenvector of AI. c) Majorana den¬ 
sity function in a finite nanowire with = 20Lu calculated 
by direct diagonalization of the Hamiltonian using covari¬ 
ant derivative discretization. The three panels correspond 
to As = 10 Eu, 0 = 75° and the rest of parameters as in Fig. 

S:. 

A. Decay lengths 

Within our complex-band-structure approach to the 
semi-infinite nanowire we can estimate the length of the 
Majorana decay tail from the imaginary part of the al¬ 
lowed wave numbers. The lower the imaginary part, the 
longer the Majorana decay tail (and thus the required 
length of the nanowire to contain it without distortion). 
Figure S] shows a typical evolution of the wave numbers 
(red islands) in the complex plane as the polar angle is ap¬ 
proaching the critical value. In the sequence from upper 
to lower panels, one of the wave numbers moves along the 
imaginary axis towards the origin; the phase transition 
being signaled by one mode touching the origin (lower 
panel). 

We calculate the required nanowire length with the 
smallest imaginary wave number of the set of all allowed 
wave numbers However, as shown in Fig. 2] the 

smallest imaginary part S(fc*'’”^) = changes from an 
approximately fixed mode to the one touching the origin 
when approaching the phase transition. We define the 
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FIG. 4. measure in a color (gray) scale showing the 

position of the allowed wave numbers as zeroes (red islands). 
The parameters used are g = 15, m* = 0.033, a = tvEuLu, 
As = 3Eu, Ab = 10 Eu- Panels from top to bottom are for 
different polar angles: 6 = 68° (a), 6 = 67° (b), 6 = 66° (c). 

mode length Lm as two times the length needed for the 
wave function to drop to one percent of its maximum, 
that is q ^n estimate of the nanowire 

length for undistorted Majoranas is simply the maximum 
of all mode lengths. 

In Fig. Owe show the mode lengths of the two allowed 
wave numbers of Fig. O As we decrease the polar angle 
from 90 degrees, the needed nanowire length (the higher 
of the two curves) remains more or less stable until 0 
approaches the critical value. A few degrees before the 
transition the Majorana contracts before diverging to in¬ 
finity at the phase transition angle. 


VI. FINITE NANOWIRES 

The phase diagram can also be checked with full di- 
agonalizations of nanowires with large, but finite, L^. 
Though more realistic, this approach is conceptually 
more qualitative, since finite-nanowire Majoranas are not 
exact zero modes but small energy modes (the smaller the 
energy the larger L^)- Equivalently, the phase bound- 



70 80 90 

0 ( degrees ) 

FIG. 5. Mode lengths Lm (defined in Sec. IV A|l for the 
two wave numbers shown in Fig. 0] Note that the required 
nanowire length at each 6 is the higher of both curves. 

aries become blurred due to the finite size effect. Figure 
[3t shows the density of the finite nanowire Majorana cor¬ 
responding to the semi infinite one of Fig.[3lD. Differences 
are small, just a slight distortion and a somewhat longer 
decay tail of the finite-nanowire density. 

In the finite nanowire diagonalization with orbital 
terms we have found it crucial using a covariant grid 
discretizatioui^Si^ Otherwise, numerical artificial biases 
wrongly suggest that Majoranas are always destroyed by 
orbital terms^^SS in clear contradiction with the phase di¬ 
agram (Fig. [2]) and the semi-infinite wire analysis. In 
essence, the covariant discretization amounts to express¬ 
ing the canonical momentum components as symmetry¬ 
like transformations. For instance. 

Although these two representations of the canonical op¬ 
erator are equivalent in the continuous limit, they are not 
on a discrete grid. 

As demonstrated in Ref. [Sll the covariant derivative 
preserves by construction the gauge invariance of the so¬ 
lutions while a non-covariant treatment only does that 
for extremely fine discretizations, unfeasible in our case. 
Changing the gauge origin usually constitutes a severe 
difficulty for numerical discretizations not using covariant 
derivative formulations. In our case, we can introduce an 
arbitrary gauge center yc for the canonical momentum, 
generalizing Eq. (ITOl) to 

^ (20) 

where the isospin sign s,- = ± is introduced in order to 
preserve the particle-hole symmetry of the Bogoliubov- 
deGennes equation. We have checked that our numerical 
results for the finite nanowire diagonalization as, e.g., in 
the lower panel of Fig. [31 do not depend on the choice 
of j/c, thus proving the gauge invariance of the finite sys¬ 
tem results. We have also obtained good agreement of 
the finite nanowire diagonalizations and the results of the 
semi-infinite system regarding the existence or absence of 
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a zero mode in the different regions of the phase diagrams 
(Fig. [2), again proving the reliability of the method. No¬ 
tice that the semi-infinite solution, being purely ID, can 
be obtained in very dense y grids, while the finite system 
2D diagonalization requires much coarser xy grids. 


VII. CONCLUSIONS 



In this work we have shown that the orbital motions 
caused by perpendicular components of the magnetic 
field in planar 2D nanowires give rise to a rich phase 
diagram, with regions containing Majoranas for sizable 
vertical tilting of the magnetic field. In fact, with proper 
parameters, it is possible to find Majoranas in a fully per¬ 
pendicular field. We have developed a general numerical 
method to obtain the Majorana phases in nanowires in a 
computer efficient way and we have checked this method 
against alternative calculations for semi-infinite and fi¬ 
nite nanowires. Analytical expressions of the transition 
boundaries in asymptotic regions have been found. For 
realistic parameter values (weak a) these analytical ex¬ 
pressions are a good approximation in general and not 
only asymptotically. In the strong orbital limit the criti¬ 
cal angles are independent of sample details. Finally, the 
relevance of the covariant grid discretization for the finite 
nanowire diagonalization has been pointed out. 
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FIG. 6. a) Phase transition boundaries for a 2D planar 
nanowire with a longitudinal (x) magnetic field. We have as¬ 
sumed As = 3Eu and y — 0. The shaded regions correspond 
to topological phases with a Majorana zero mode. 


interaction can be absorbed in the kinetic term 



ma^ 


-I- As Ta: -I- Ab<7x , 


(A.2) 


where Py = Py — maax/h. Using a basis of square-well 
eigenstates of energies e„ = h'^Tr'^n^/{2mLy), the matrix 
to diagonalize is 



(A.3) 


The diagonalization of this simplified Hamiltonian yields 
the eigenenergies 


Er, 


SlAg + S2 



ma^ 


+ A^ 


(A.4) 
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Appendix: Analytical limits 
1. Longitudinal magnetic field 

When the magnetic field is along x (see axis orienta¬ 
tions in Fig. 1) the phase transition law is fully analytical. 
As discussed in Sec. IIIIl finding the phase transition im¬ 
plies searching for the zero energy eigenstates of the sim¬ 
plified Hamiltonian h given in Eq. For an cc-oriented 
field = 0 and all orbital terms vanish. Assuming also 
a vanishing bottom potential for V(y) we have 

h = - ^^T^+ As Tx + (^Ab - o’x • (A.l) 


with n = 1,2,..., Si = ±1 and S 2 = ±1. Of the four 
eigenenergies only the two with opposite si and S 2 can 
lead to a zero energy solution at the critical values 

^b]u = + + A2 . (A.5) 

Equation (IA.5I) with n = 1,2,... gives the critical Zee- 
man parameter of phase transitions for a two dimensional 
nanowire in parallel magnetic field. Notice that in Eq. 
(IA.5I) n has to be interpreted simply as an ordering in¬ 
dex of the successive transitions, and not as a label of 
independent transverse modes. These latter interpreta¬ 
tion would be wrong, since different transverse modes are 
coupled through the Rashba mixing term and one can not 
associate a particular transverse mode with a particular 
transition point. Shaded regions in Fig. [5] contain one 
Majorana mode, while white regions have none. There 
are no regions with multiple Majoranas due to the en¬ 
ergy splittings induced by the Rashba mixing in planar 
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nanowires . 1 

2. Strong orbital limit 


The eigenstates of Eq. (lA.ll) can be obtained analyt¬ 
ically noticing that the linear py term from the Rashba 


When the kinetic orbital effect overcomes both the con¬ 
finement by the transverse square well and the Rashba 



























Panel 


a) 

1.43 

0.12 

0.001 

0.36 

b) 

1.03 

0.12 

0.48 

0.36 

c) 

0.74 

0.12 

1.22 

0.36 

d) 

1.06 

0.12 

1.98 

0.36 


TABLE I. Numerical values in effective units of the inequal¬ 
ities, Eq. (11) of the paper, corresponding to the four panels 
in Fig. 2 of the paper. 


term, the magnetic length = \JhcleB^ is smaller than 
ly (here we define ly = Ly of Fig.[T]) and also smaller than 
the Rashba length 1^ = T?Ima. In the limit <C (?y, la) 
it is possible to derive an analytical expression of the 
transition boundaries. Neglecting the square well V(y) 
and the Rashba terms in Eq. (jH]) we find 


h = 




2m 2mlt 


pI 


- p]tz +Asa-h +As Tj: . (A.6) 


The eigenvalues of Eq. (IA.6I) are straightforward in a 
basis IuSo-St-), where n = 1,2,... represent now harmonic 
oscillator eigenstates, = ± indicates spin eigenstates 
in direction n, while Sr = ± indicates isospin in direction 
z. Since the h matrix is diagonal in spin, we can diag¬ 
onalize each subspace independently. Eor instance, the 
matrix for So- = + reads 


(ho) 


— P + Ag 

As 


As 


(e^)-/r) +Ab j ’ 


(A.7) 


with = {n — \l2)Pi} jml^. The eigenvalues of Eq. 
(IA.7I) are easily found, as well as those of the analogous 
matrix for spin = —. 


The null-eigenvalue condition for h is now 


Ag = 





n 2 


■AL 


(A.8) 


that looks similar to Eq. (IA.5I) . An essential difference, 
however, is that the r.h.s. in Eq. (IA.8I) depends itself on 
the Zeeman parameter through It is 


ml"]. 


- Ag cos6 , 

gm* 


(A.9) 


where m* is the ratio of effective to bare mass, m = 
m*me, while g is the gyromagnetic factor defined from 
the Zeeman parameter by As = gygB/2. Erom Eq. 
(IA.8I) we finally arrive at the following relation 


4 („-l)A 


Eor large enough Ag, as compared to Ag and y,, this leads 
to the prediction of field-independent critical angles 


cos 




gm 


(A.ll) 


_ 4(n- i) ’ 

as given in Eq. (Il4ll . 

The triple inequality Iz {ly,la,ls), where we define 
Is = /mAg, leads, when written in effective units, 
to Eq. (na. In this situation the phase diagram does 
not deviate much from the straight lines of the analytical 
limits, Eqs. (IA.5I) and (lA.llI) . Table U contains the nu¬ 
merical values of the inequality sides for the four panels 
in Fig. 2. While panel a) fulfills all conditions, for the rest 
of panels the second inequality degrades as a increases 
from panel b) to d). This explains the deviations in those 
panels from the analytical limits. 
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